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We analyze the dynamics of a single scalar field in Friedmann-Robertson- Walker universes with 
spatial curvature. We obtain the fixed point solutions which are shown to be late time attractors. 
In particular, we determine the corresponding scalar field potentials which correspond to these 
stable solutions. The analysis is quite general and incorporates expanding and contracting universes 
with both positive and negative scalar potentials. We demonstrate that the known power law, 
exponential, and de-Sitter solutions are certain limits of our general set of solutions. 
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I. INTRODUCTION 



Scalar fields have played a very important role in mod- 
ern cosmology. Today's observed acceleration of the uni- 
verse, for example, may be explained by the dynamics 
of a scalar field (for a review, see [H). The scenarios 
proposed to solve the initial conditions of the standard 
Big Bang theory such as inflation [3;^ H , pre-big bang 
[liMIll, and the ekpyrotic/cyclic 0,13j1^ scenarios also 
usually require a scalar field. What makes these scenar- 
ios interesting is the existence of attractor solutions which 
implies that the dynamical system becomes insensitive to 
the initial conditions. A class of such attractor solutions 
is referred to as scaling solutions, where the energy den- 
sity can be divided into contributions which scale with 
one another throughout the evolution of the universe. In- 
vestigating the nature of scaling solutions allows one to 
understand the asymptotic behaviour of a particular cos- 
mology and helps determine whether such behaviour is 
stable or just a transient feature. 

In a spatially flat Friedmann-Robertson- Walker 
(FRW) universe filled with a perfect fluid and a canon- 
ical scalar field, such scaling solutions are obtained 
through a simple exponential potential The at- 

tractor behaviour in this system has been analyzed ex- 
tensively [m, [H, [H, [13, [ll, More recently they 
have been obtained for a wide class of modified cos- 
mologies proposed by fundamental theories. The cos- 
mologies inlcude those with a non-minimally coupled 
scalar field [TtI. [!§. [TqI. [20|. [21|. [2^ . braneworld cos- 
mologies tachyon cosmologies 



I, tachyon cosmologies 

, ii [MJliS m, [H, [3, loop 

smology, 4ll |42| . l43l[rp hantom cosmology [ij 



quantum cosmolog 

and Gauss-Bonnet cosmology [45|, [46| (For examples of 
similar analysis in the context of modified gravity, see 
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[£71, [48|, [4g, [50, [Ml, [52, [Sa [5J, |55| ) . 

Fundamental theories may introduce modifications to 
the Friedmann equation in various models that may or 
may not be of a similar form. Since the stability of 
the background cosmology is understood to be a desir- 
able property, the pragmatic approach is to accommo- 
date such modifications in the phase space analysis of 
the standard dynamical system. For the case of a flat 
FRW universe, the analytical method of obtaining the 
scaling solutions has been extended to a general set-up, 
where the Hubble parameter is given either by an arbi- 
trary power [sgI [57[ or an arbitrary function 58] of the 
total energy density. 

According to the latest results of WMAP5 al- 
though limits being placed on the flatness of our ob- 
servable universe are sqeezing the possible outcomes, 
-0.063 < ilk < 0.017, there is stfll the possibility that 
we live in a spatially curved FRW universe. It is ap- 
propriate, therefore, to investigate the consequences of 
curvature on the stability of the background cosmology. 

In this paper, using the method developed in [s^ in the 
context of braneworld corrections made to the standard 
cosmology, we will obtain the scaling solutions in spa- 
tially curved FRW universes and classify the asymptotic 
behaviour of the systems. We aim to present a general 
approach to analysing such systems by deriving the form 
of the scalar potential leading to late time attractors. 
In doing so, we will demonstrate that our results reduce 
to the known results in the literature once applied to 
the corresponding models [13, [6l[- (considerable atten- 
tion has also been given to scaling solutions in Bianchi 
spacetimes [H, [H, [H, [H, [H, [131). As was done in 
for the case of a flat FRW universe, we consider positive 
and negative potentials of the scalar fleld, since negative 
potentials can provide interesting cosmological scenarios 
both in an expanding 143,, ^1 and collapsing universe 

Big. 

The rest of the paper is arranged as follows. In Sec. II 
we present the equations of motion and introduce the 
variables which allow the scaling solutions to be deter- 
mined. After analyzing the stability of these solutions 
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in Sec. Ill, we obtain the general relations which hold 
when the scaling solutions are obtained in Sec. IV. Then, 
we apply these results to the open and the closed FRW 
universe in Sec. V and Sec. VI, respectively. Finally, we 
summarise in Sec. VII. 



II. EQUATIONS OF MOTION 

We consider Friedmann- Robertson- Walker (FRW) cos- 
mologies such that the dynamics is determined by an ef- 
fective Friedmann equation of the form 



(1) 



where H = a/a is the Hubble parameter, a is the scale 
factor, p{a) is the total energy density of the universe, 
a dot denotes differentiation with respect to cosmic time 
and TO4 is the four-dimensional Planck mass. Through- 
out this paper, we work in units where 87r/m| = 1, and 
we introduce the notation H = esfW^, where e = ±1 
corresponds to an expanding or a contracting universe, 
respectively. Modifications to standard relativistic cos- 
mology are parameterized by the function L{p{a)) and 
this is assumed to be positive-definite without loss of 
generality. 

We will investigate models where the universe is 
sourced by a self-interacting scalar field cf) with poten- 
tial V{<j)) together with a barotropic fluid with equation 
of state p.^ = (7 — where 7 is the adiabatic index. 

The energy density and pressure of the scalar field are 
given by — 0^/2 -I- V and — ^ jl — V, respectively. 
We note here that the effective adiabatic index of the 
scalar field is given by 7^ = (p^ +p^)/pcf,. As in conven- 
tional cosmologies, we assume that the energy-momenta 
of these matter fields are covariantly conserved and this 
implies that 



(5) 
(6) 



Yn = -e\^^XY+^Y[2X^+j{l-X^-aY^)] 



X j^ = -eV6X'^{T-l)X + 3X[2X'^+j{l-X^ -aY"^)] 
d\uL 



P- 



dp ' 



(7) 



where A'' = Ina, the subscript N denotes differentiation 
with respect to this parameter, and 



(8) 



For this new set of variables, the definition of the total 
energy density implies the constraint equation 



X^ + aY^ + ^ = l. 



(9) 



From Eq. ([2]) the scale factor, a, is expressed in terms 
of the fluid energy density, p^, as 



1 1^ 



(10) 



where the quantities with subscript («) are evaluated at 
some initial time. 

Differentiating A, given by Eq. with respect to the 
scalar field, (j), for the special case of A being a constant, 
one can show that 

„ dlnL , ^ 

This is the case we consider for the rest of the paper. 



p^ -SjHpy, (2) 
4> = -3H^-V^, (3) 

where the subscript (f> means differentiation with respect 
to the field. Eqs. (HJ-® close the system that deter- 
mines the cosmic dynamics. Introducing the variables 

[mi 



III. STABILITY 

In order to carry out the stability analysis for this case, 
it is sufficient to solve Eqs. (O-® for variables X and 
Y, since Eq. ([7]) reduces to the constraint equation (|lip . 
For this system of equations we find the following set of 
physical fixed points 



X 



2p 



Y 



(4) 



where p is the total energy density of the universe, and 
we adopt the notation F = for a = ±1 for positive 
and negative potentials, respectively. Eqs. (P)-© can be 
rewritten in the form 




(12) 
(13) 

(14) 
(15) 
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In what follows we describe the classification of these 
solutions in terms of their region of existence and stability 
for an expanding universe containing a scalar field with a 
positive potential (i.e. e = 1, and a = 1). Later, we will 
generalise our analysis for an expanding or a contracting 
universe in which the scalar field could have cither a pos- 
itive or a negative potential. These are summarised in 
Table D 

The first two points in (|12p correspond to the scalar 
field kinetic energy dominated solutions. The third point 
corresponds to the fluid dominated solution. Upon 
imposing the constraint Eq. ([5]), the fourth point is 
a solution that exists for A 0, < 7 < 2, and > 
37; and it describes a scenario where, for a given fluid, 
the contribution of the scalar field density to the total 
energy density scales with that of the fluid to the total 
density, i.e. Xl + aY'^ — For convenience, we refer 
to this solution as the fluid-scalar field scaling solution 
throughout this paper. The final point arises if A^ < 
6, and the constraint Eq. ^ implies that, in this case, 
the energy density of the universe is dominated by the 
scalar field, having an effective adiabatic index jtt> — 
This describes a scaling solution, where as the universe 
evolves, the kinetic energy and the potential energy of 
the scalar field scale together. We refer to this solution 
as the scalar field dominated scaling solution throughout 
this paper. 

Having obtained the scaling solutions, we need to in- 
vestigate their stability to small fluctuations. Consider- 
ing perturbations of the form 

X ^Xc + 6X, Y = Yc + SY , (16) 

where SX oc e""^, and SY oc e'^^ , and expanding Eqs. 
(in])-®, for the kinetic energy dominated solutions, 




respectively. In this case, stability is only achieved for 



fluids with adiabatic index of 7 > 2, which is not satis- 
fied for known realistic fluids. These points, are therefore 
considered to be unstable. For the fluid dominated sce- 
nario, the eigenvalues are given by 



which is clearly only stable for negative values of the adi- 
abatic index, 7. We will therefore consider this solution 
to be unstable for all realistic types of fluid. The fluid- 
scalar field scaling solution yields the eigenvalues (recall 
< 7 < 2, and A^ > 87) 

which implies unconditional stability when these solu- 
tions exist. The same analysis yields the following eigen- 
values for the scalar field dominated scaling solutions 



w+ = -{X^-6) w.=X^-3j, (20) 

which indicates that when these solutions exist, as long 
as A^ < 87, stability is guaranteed. We also note that 
in a contracting universe (e = —1), an attractor solution 
correpsonds to one with positive real eigenvalues. This 
is because the parameter with respect to which the dy- 
namical system of (H])-® is described, N, becomes a de- 
creasing function of time. Keeping this in mind, the gen- 
eral consideration of various combinations of an expand- 
ing/contracting universe containing a scalar field with a 
positive/negative potential are captured in Table U be- 
low. For the case of L = 1 (i.e. a flat FRW universe), 
a similar classification was done for simple exponential 
potentials in p^ . 



J 





X = 1 
y = 


X ^-1 
y = 


X = 
y = 


y = v/^a 






exists 


stable 


exists 


stable 


exists 


stable 


exists 


stable 


exists 


stable 


e = +1 ,a = +l 


VA, V7 


No 


VA, V7 


No 


VA, V7 


No 


A / 0, A" > 37 


when exists 


A" < 6 


A" < 37 


6 = 4-1 , a = -1 


VA, V7 


No 


VA, V7 


No 


VA, V7 


No 


No 




A" > 6 


No 


e = -1 , a = +1 


VA, V7 


> A > -Ve 


VA, V7 


o<x<Vg 


VA, V7 


No 


A / 0, A" > 37 


No 


A" < 6 


No 


e = —1 , a = —1 


VA, V7 


> A > -Vg 


VA, V7 


o<\<Vg 


VA, V7 


No 


No 




A" > 6 


when exists 



TABLE I: This table summarises the existence and stability conditions for an expanding (e = 1) or a contracting (e = — 1) 
universe containing a fluid with the adiabatic index 7 and a scalar field with either a positive (a = 1) or a negative (a = — 1) 
potential. 



IV. GENERAL RELATIONS FOR SCALING scalar field potentials which correspond to these late time 

SOLUTIONS 



In the remainder of this paper we concentrate on the 
solutions given by Eqs. (fT4| -(|15 p . We aim to derive the 
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attractors for particular forms of L{p) based on spatially 
curved cosmologies. Working at the fixed points, from 

Yc — y^^, and the fact that = X"^ + aY^ = con- 
stant for these solutions, L{p) can be described as L{V). 
Moreover, from the definition of A, one can show that 



A. Case A: Fluid-scalar field scaling solutions 

When the universe is expanding (e — 1), the solution 
given by ([14]) exists and is stable for > 87, and the 
scale factor a can be expressed in this case as 



jwL-'*^ (21) 

Given a specific form of L, therefore, it is possible to 
derive the scalar potential resulting in the scaling solution 
by integrating Eq. (PTjl . We note that, as demonstrated 
in Eq. (j2ip is equivalent to imposing the constraint 
Eq. (fTTj) throughout the evolution of the field. 

Furthermore, since at the fixed points, X = Xc, (j) is 
a monotonically varying function of cosmic time, t, and 
can be considered as a suitable dynamical variable for 
the system. We note here that this assumption is only 
valid when Xc 7^ 0. In general, the scalar field Eq. ([3]) 
can be expressed in the form 



a = Ap-^^ , with A = a(,)p4) (a^) ' ^^S) 

The correction function L given by Eq. ([25]) can then be 
rewritten as 




By considering the fixed point, substituting this into 
Eq. ([21]) and integrating, yields the scaling solution po- 
tential 



P0 = -mi\? , (22) 

and for the scaling solutions, this equation can be ex- 
pressed as 

but since V and (/) are related through Eq. (|2ip. this can 
be integrated to find i as a function of the scalar field 

t^~aeY,^fl{^ f d(\>L{V)^^. (24) 

V. OPEN FRW UNIVERSE 

Here, we consider what form of the scalar field po- 
tential provides the fixed point solution characterised by 
{X,Y) — {Xc,Yc) in the open FRW universe. In this 
case, L{p) is given by 

L(p) = ./l + A. (25) 

In our analysis, we assume that A ^ 0, but we return 
to this point and consider the special case of A = at the 
end of this section. We note at this point that an expand- 
ing (contracting) open universe obeying the Friedmann 
equation ([l} could stop its expansion (contraction) pro- 
cess and begin to contract (expand) if the scalar field has 
a negative (positive) potential. It is, therefore, interest- 
ing to study the physics obtained from different combi- 
nations of e and a in an open universe. We begin our 
discussion by concentrating first on an expanding (e ~ 1) 
universe sourced by a positive potential (a = 1) scalar 
field. 



|y(0)r = Y^>^ j cosech^ [^-pct^^ , (28) 

where /i = (2 — 37)737. Notice that the valid range of p, 
for the region of existence of these solutions, is —2/3 < 
/i < and < p < 00. The special case oi p — 0, where 
Eq. ([^5]) is no longer valid, is discussed later. 

Using the definition of Y at the fixed point, and substi- 
tuting Eq. ([28]) into Eq. ([27]) , the correction function can 
be written in terms of the scalar field as L(0). Upon sub- 
stituting this form back in Eq. ([M]) the time dependence 
of the scalar field can be evaluated as 

Classifying the behaviour of the evolution equations 
above in terms of the sign of the parameter p^ and recog- 
nising the sign of the argument inside the brackets re- 
mains invariant under the transformation A — > — A and 
(f> —(j), we choose to work in the first quadrant without 
loss of generality. 

For /i > 0, we find that at early times, as (/) — > 0, the 

asymptotic form of the potential relating to the curvature 

— 2. 

dominated universe is a power law function V (yi (j) 1^ . 
However, at late times, as (f) 00, the potential of a 
fluid-scalar field dominated universe is of an exponential 
form V ^ exp[— A(/)]. 

On the other hand, for p < 0, at early times, as 
— > cx), where the curvature is negligible, the asymp- 
totic form of the potential is obtained to be an expo- 
nential one V exp[— A(/)]. Once the universe becomes 
dominated by the curvature at late times, as (/> — > 0, the 
potential can be approximated by a power law function 
V (X (f) , recovering the standard result a ^ t. 
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For /X = 0, since the correction function L is a constant 
and can be thought of as modification to Newton's grav- 
itational constant in Eq. ([1]), the expansion law becomes 
that of the flat Friedmann cosmology, and the potential 
yielding this late time attractor will then be the exponen- 
tial potential as found in [lB| • We further note that 
this scenario refers to the special case where the contri- 
bution of the scalar field, the fiuid, and the curvature to 
the Friedmann equation scale together. 

We note, according to the results illustrated in TablelH 
that the fixed points do not exist for negative poten- 
tials. This is true for an expanding or a contracting uni- 
verse. Although in an expanding (contracting ) universe, 
a negative (positive) scalar potential could slow down the 
growth of the scale factor and cause the universe to col- 
lapse (expand), we do not expect the late time attractors 
to be given by these scaling solutions. We can also see 
from Table [T] that this set of solutions does not exist in 
a contracting universe sourced by a negative scalar field 
potential. 

B. Case B: Scalar field dominated scaling solution 

When the universe is expanding (e = 1) the solution 

given by p5|) . is an attractor for < 87, and the scalar 

field dominates over the fluid, resulting in the effective 
\ 2 

adiabatic index 7^ = The scale factor is then given 
by 

a = Bp-* , with B = a(,)P(^ . (30) 



The form of the correction function L given by Eq. 
([25| can then be rewritten as 



L{p)=^l + ^p'i^. (31) 

By considering the attractor solution, substituting this 
form of the correction function into Eq. (HI]), and inte- 
grating, yeilds the scaling solution potential 



where = (2 — A^)/A^, and the form of the potential is 
valid for A V 2 . 

By using the definition of Y at the fixed point, and 
substituting Eq. ([32]) into Eq. ([3T|) the correction function 
can be written in terms of the scalar field as L{(f)). From 
this form the time dependence of the scalar field can be 
evaluated using Eq. ([24l) to be 

^ = (|,)^/^0sinhV^(^0), (33) 

By analogy, the asymptotic behaviour of the universe 
in this scenario can be obtained using the same method 
as employed in Case A. These results are summarised in 
Table n 





Case A 


Case B 




M < 


^ = 


M > 


u <Q 


u = 


> 


Early times 


V ~ exp[-A<?!>] 

negligible 
curvature 


V ~ exp[-2(j!>] 

scaling 
curvature 


V ~ (^~~ 

curvature 
dominated 


V ~ exp[-A(^] 

negligible 
curvature 


V ~ exp[-2(^] 

scaling 
curvature 


V ~ <l)-i 

curvature 
dominated 


Late times 


V ~ (p'~ 

curvature 
dominated 


V ~ exp[-2(/>] 

scaling 
curvature 


V ~ exp[-A<?!>] 

negligible 
curvature 


V ~ (j>-^ 

curvature 
dominated 


V ~ exp[-2(^] 

scaling 
curvature 


V ~ exp[-A(^] 

negligible 
curvature 



TABLE IL This table summarises the asymptotic behaviour of an expanding (e = 1) open FRW universe described by the 
scaling solutions, when the scalar potential is positive. Case A refers to the fluid-scalar field scaling solution, and Case B 
corresponds to the scalar field dominated scaling solution. 



We now turn to our general set of results in Table U 
and note that in a contracting universe (e — —1), if the 
scalar potential is negative {a = —1) and steep enough 
(A^ > 6), the scalar field dominated scaling solution can 
also result in stable attractor behaviour. Such solutions 
will be unconditionally stable if they exist, and we find 



the form of the potential together with the time evolution 
of the scalar field are still given by Eqs. ([32]) and ([33| . In 
this scenario, u is clearly negative, and due to the time 
reversal of Eq. for e = — 1, one expects from Table HIl 
an early time power law behaviour of the potential to be 
followed by an exponential form at late times. This is the 
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generalisation of [15| to general curved space scenarios. 

C. Case C: A ^ 

We now return to the special case of A = 0, and no- 
tice that this corresponds to a constant potential (i.e. 
de-Sitter space). Solutions for this type of universe are 
known (70| , and it is therefore interesting to see how they 
would fit into the larger class of solutions we are present- 
ing here. Considering the scaling solutions, we notice 
that the fixed points described by (|14p do not exist for 
the class of constant potentials. However, the attrac- 
tor solutions (fT5|) reduce to = 0, and Yc = 1, which 
describe an exact de sitter solution. As mentioned previ- 
ously, the argument of using </) as a monotonically varying 
function of time breaks down when Xc oc (j) — 0. For this 
simplified case, one can clearly solve Eqs. (IlJ-© directly 
to find the exact solution. We notice, however, that we 
should also be able to recover this solution by investigat- 
ing the behaviour of the scalar potential and the scale 
factor very close to the fixed point. We do this by con- 
sidering small values of A. From Eq. ([55]) . we find the 
asymptotic dependence of the scalar field to be given by 
t oc (j). Substituting this into Eq. ([32|) . and using Eq. 
(fTO|) near the fixed point, we find that the scale factor 
evolves as a oc sinh (-^i^i). Thus, recovering the de-sitter 
solution as a special case of our larger set of solutions 
derived here. 

We further note from Table U that when A w 0, a neg- 
ative potential scaling solution docs not exist in an ex- 
panding or a contracting universe; and the positive po- 
tential solution, which exists in a contracting universe, is 
unstable. 

VI. CLOSED FRW UNIVERSE 

Here, we aim to derive the forms of the scalar poten- 
tials which would result in stable fixed points p^ - p^ . 
corresponding to the scaling solutions in a closed FRW 
universe. In this case, L{p{a)) is given by 

L{p) = Jl-^. (34) 

In this scenario, the curvature can not dominate the 
contributions to the total energy density of the universe, 
otherwise the correction function ()34|) would be imagi- 
nary and the right hand side of the Friedmann equation 
would become negative. The valid range of the correc- 
tion function is thus given by < L < 1. The case of 
L — Q corresponds to a universe with a constant scale 
factor, and describes the case where there is no expan- 
sion or contraction taking place. The universe may un- 
dergo a bounce at this point and change from an expand- 
ing (contracting) behaviour to a contracting (expanding) 
one. We bear this in mind in our following analysis. As 



in the open case, we first concentrate on the expanding 
universe (e = 1) with positive potential (a = 1) scalar 
field, and we will comment on the negative potentials and 
the difference between an expanding and a contracting 
universe for each set of scaling solutions. 

A. Case A: Fluid-scalar field scaling solution 

In an expanding universe (e = 1), the solution given by 
is an attractor for A^ > 37, and describes a scenario 
where the contribution of the scalar field scales with that 
of the fluid. The scale factor in this case is also given by 
Eq. ((26|) . and after substituting this into Eq. (|34|) . the 
form of the correction function, L, is found to be 

Notice there is a maximum level of density, that de- 
pends on A, beyond which L lies outside its valid range 
and the situation is unphysical. Following our previ- 
ous analysis, except for the case p, — Q, where /i = 
(2 — 37)737, the potential yielding the fixed point so- 
lution is given by 

in0)r =i;'^(^)sech2QM0) • (36) 

Using Yc — \/W\/P: ^iid substituting Eq. ([36]) into 
Eq. ([35|) . the correction function can be written as L{(j)) = 
tanh . And substituting this form of L back into 

Eq. 1^4)) . yields the time dependence of the scalar field 
as 

In order to examine the asymptotic behaviour of the 
scalar field potential, we follow our previous line of ar- 
gument and divide up the region of validity of /i into its 
positive and negative values ( i.e. — 2/3 < /i < 0, and 
< fj, < 00). Once again, making use of the symme- 
try A — A and — > — 0, we restrict our analysis to 
the potentials lying in the first quadrant, without loss of 
generality. 

For /i > (i.e. 37 < 2), in an expanding universe 
at early times, since > a'^'^, one would expect the 
curvature to dominate initially. However, as mentioned 
above, this can not happen in a closed universe, oth- 
erwise the correction function (j35p becomes imaginary 
and consequently, the right hand side of the Friedmann 
Eq. ([T|) will be negative. A pragmatic starting point for 
the expanding evolution is, therefore, where the curva- 
ture contribution is only just subdominant to the scalar 
field and the fluid densities. We set this point to cor- 
respond to t « 0, since going back in time from this 
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point, in the absence of motivations from a fundamen- 
tal theory, is not physically meaningful. This is when 
the total energy density approaches its maximum value 
P ~* Pmax = and so, the correction function 

L{(p) — > 0. Eq. Il]) then suggests that the scale factor 
does not change as ^ 0. Once the universe starts 
expanding, the contribution of the fluid and the scalar 
field to the energy density starts to dominate that of the 
curvature. In the limit where the curvature has become 
negligible and the universe asymptotes to a flat FRW 
spacetime, we recover the exponential limit of the poten- 
tial V ~ exp[— A(/)]. 

For /J, < 0, as the universe expands, it starts from an 
asymptotically flat FRW state dominated by the fluid 
and the scalar field. The potential in this case is that of 
the exponential form V ^ exp[— A(/)], and the curvature 
can be neglected. In time, we reach the maximum level 
of energy density, where the contribution of the curva- 
ture becomes very close to that of the combination of the 
scalar field and the fluid, and the scale factor seizes to 
grow as H ^ 0. After this turning point, the universe 
begins to collapse, and this provides us with a natural 
motivation to consider a contracting universe. However, 
one can see from Table HI that in a contracting universe, 
late time fluid-scalar field scaling attractor solutions ei- 
ther do not exist at all (for a = — 1) or are uncondition- 
ally unstable when they do exist (for a = 1). Once the 
universe starts to collapse, solutions will asymptote to- 
wards the kinetic dominated solutions for < 6. On the 
other hand, for A^ > 6, solutions will asymptote towards 
the fluid dominated scaling solutions if the potential is 
negative {a = —1), and will be unstable for positive po- 
tentials. 

It is worth noting that in the limit of H — 0, as a 
result of L 0, the total energy density is a constant 
close to its maximum value. From Eqs. dH) and ((36)l . 
we conlcude that in this limit the potential is almost 

constant V ^ j . One may then naively expect 

the universe to expand exponentially as it would do in a 
de-sitter case; however, we note that the correction factor 
stops this behaviour by keeping the scale factor constant 
in Eq. (P). 

The special case of ^ = corresponds to where the con- 



tribution of the curvature, the fluid, and the scalar field 
are all scaling at all times during the evolution. We note 
that this is only valid if the curvature does not dominate 
over the combination of the scalar field and the fluid, in 
which case, the correction function of Eq. ([35]) is a con- 
stant and the potential has the exponential form. 

Similarly to the open universe scenario described ear- 
lier, we can see from Table U] that for negative scalar po- 
tentials {a — —1) this set of scaling solutions do not exist 
in either an expanding universe or a contracting one. 
B. Case B: Scalar field dominated scaling solution 

In an expanding universe (e = 1) when the potential of 
the scalar field is positive {a = 1), the fixed points given 
by are attractor solutions for A^ < 37, and the scale 
factor evolves according to Eq. ([30| . Therefore, the form 
of the correction function given by Eq. (|34p becomes 



Lip) = ^l-^p'^. (38) 

We remember that in this case, the contribution of the 
fluid is negligible for these solutions at all times during 
the evolution of the universe. 

Similarly to our previous analysis, we find the corre- 
sponding scalar potential, except for the case of v — 0, 
to be given by 

in</')r = >;'''(|^)sech2Qz.</.) , (39) 

and the time dependence of the scalar field is 

By analogy, the asymptotic behaviour of an expanding 
universe can be obtained through a similar discussion as 
the previous case. These results are summarised in Table 

imi 



It is worth noting at this stage that according to the 
results demonstrated in Table HI an expanding (contract- 
ing) universe with a negative (positive) scalar potential 
does not admit attractor solutions. However, for a con- 
tracting universe with a negative potential, if such solu- 
tions exist they will be unconditionally stable, and the 
form of the potential together with the time evolution of 
the scalar field can be seen to be given by Eqs. (|39p and 



601). 



C. Case C: A ^ 

As mentioned in the open universe scenario, this case 
yields an ill-defined solution corresponding to the fixed 
points described by (|14p . However, A = results in a 
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Case A 


Case B 




M < 




M > 


V <Q 






Early times 


V ~ exp[-A(^] 

negligible 
curvature 


V ~ exp[-2<^] 

scaling 
curvature 


curvature 
only just 
subdominant 


V ~ exp[-A<^] 

negligible 
curvature 


V ~ exp[-2<?!)] 

scaling 
curvature 


curvature 
only just 
subdominant 


Late times 


v-=»i-.' (*)"■■ 

curvature 
only just 
subdominant 


V ~ exp[-2<^] 

scaling 
curvature 


V ~ exp[-X(t>] 

negligible 
curvature 


y.r^(f)"' 

curvature 
only just 
subdominant 


V ~ exp[-2,^] 

scaling 
curvature 


V ~ exp[-A-^] 

negligible 
curvature 



TABLE IIL This table summarises the asymptotic behaviour of an expanding (e = 1) closed FRW universe described by the 
scaling solutions, when the scalar potential is positive. Case A refers to the fluid-scalar field scaling solution, and Case B 
corresponds to the scalar field dominated scaling solution. 



perfectly well-defined scaling solution when substituted 
in the fixed points given by (jlSp . The scaling solution 
then reduces to Xj. = 0, and Yc = 1, which is clearly a 
de-Sitter solution. Once again, the assumption of <f) being 
a monotonically varying function of time breaks down 
at Xc = 0. Thus, we take the pragmatic approach of 
investigating the form of the potential and the behaviour 
of the scale factor near the scaling solution by considering 
small values of A. We then find the potential to be given 
by y (X cosh^^^"^ (f'^'^) ' ^'^'^ time dependence of the 
scalar field to be t cx which from Eq. ([30]) . yields the 

scale factor as a (x cosh {^^vt^ ■ This is the form of the 
scale factor evolution in a closed de-sitter universe, and 
is a special subset of the solutions we have derived here. 

In an analogous way to our description for an open 
scenario, the fixed points IjlSp do not exist in an anti- 
de-Sitter universe; and even though they exist in a con- 
tracting universe sourced by a scalar field with a positive 
potential, such solutions are unconditionally unstable. 

VII. SUMMARY 

In order to understand the asymptotic behaviour of 
a particular cosmology and determine whether such a 
background is stable or not, a class of attractor solutions 
referred to as the scaling solutions play an important 
role. By now, even though many scaling solutions have 
been obtained in various cosmologies, they are limited to 
the case of a flat Friedmann-Robertson- Walker (FRW) 
universe. In this paper, we have analyzed the dynam- 
ics of a single scalar field in FRW universes with spatial 
curvature. We started by generalising the approach de- 
veloped in [ssl to incorporate expanding and contracting 
universes filled with a perfect fluid and a single scalar 
field. Due to a growing interest in negative scalar po- 
tentials and the possibility of obtaining them from fun- 
damental theories, these are also accommodated in our 
formalism. 

We have identified two sets of scaling solutions. One, 
where the contribution of the fluid to the total energy 



density scales with that of the scalar field throughout the 
evolution, and the other, where the scalar field dominates 
over the fluid, and the kinetic energy of the field scales 
with its potential energy. By concentrating on both types 
of scaling solutions, we obtained the generalised dynami- 
cal equations. Once a particular form of the modification 
function is given and its dependence on the total energy 
density is known, these equations can be used to derive 
the form of the scalar potential that leads to the late time 
attractor solutions. This analysis is explicitly carried out 
for the cases of an open and a closed universe. 

After presenting the general form of the potential, we 
examined the asymptotic behaviour of the dynamics. In 
an open universe sourced by a positive scalar potential, 
we concluded that in regions where the curvature is dom- 
inant, the potential can be approximated by a power law 
function; and where it is negligible, an exponential form 
is a good approximation. This result is consistent with 
[6l|, where the exponential form of the potential was 
found to be unstable when the curvature term becomes 
important. We have also shown how the well known de- 
sitter solutions can be found as certain limits of our gen- 
eral potential. We then highlighted the case of a collaps- 
ing open universe (which could happen in the presence of 
a negative potential [6^) and concluded that, if the po- 
tential is steep enough, the scalar field dominated scaling 
solution is a late time attractor. This is in line with the 
known results in the context of an ekpyrotic collapse in 
the flat FRW universe UMM- 

The general form of the scalar potential is also pre- 
sented for a closed universe. Here we are facing an in- 
teresting scenario and care needs to be taken since the 
curvature is forbidden to dominate over the combination 
of fluid and the scalar field in an expanding universe when 
the potential is positive. In the limit where the curva- 
ture is comparable to this combination, the potential is 
a constant and so is the scale factor. When the curva- 
ture becomes subdominant, the potential asymptotes to 
an exponential form. In a contracting closed universe we 
find that a steep negative potential provides the late time 
attractor in the form of a scalar field dominated solution. 
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Negative exponential potentials have been known to be 
stable solutions in ekpyrotic models, and this is a special 
limit of our derived potential. 

In this work we have concentrated on the scaling so- 
lutions, but it is also worth noting that in a contracting 
universe, if the potential is flat enough (A^ < 6), the 
kinetic dominated solutions are shown to be late time 
attractors for both the positive and the negative poten- 
tial scalar fields. This solution corresponds to the ones 
obtained in the pre-big bang cosmology 1, H, [1] ■ 

Although the method proposed in [58 1 had been mo- 
tivated by generalising a large class of modifications in- 
troduced to the Friedmann equation in the context of 
braneworld cosmologies in flat space, we have demon- 
strated here that a spatial curvature term can also be 
written in a similar fashion. We recognise that at the 
background level it is not possible to distinguish between 
corrections originating from fundamental theories that 
can be encoded within the correction function, and 
information about the spatial curvature which could be 
represented through a similar modification function. In 



order to differentiate between the two, one needs to ex- 
amine signatures such as the behaviour of perturbations 
in the universe. This analysis is beyond the scope of our 
work in this paper. Finally, in this paper, we have lim- 
ited our analysis for the canonical scalar fields. We do 
not envisage the obvious extension to the non-standard 
scalar fields like non- minimally coupled fields Ts'l , tachy- 
onic fields ^31j, or phantom fields 441 be complicated. 
It may be interesting to consider the dynamics of such 
scalar fields as well as multi-field models [7lj in curved 
FRW universes. 
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